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By A. J. Eggers, Jr., David H. Dennis, 
and Meyer M. Resnikoff 


SUMMARY 


V- 

Approximate. shapes -'6f nonlifting todies having minimum pressure 
foredrag at high supersonic* airspeeds are calculated. With the aid of 
Newton* s law of resistance, the investigation is carried out for various 
combinations of the conditions of given body length, base diameter, sur- 
face area, and volume. In general it is found ' that when body length is 
fixed, the body has a blunt nose; whereas, when the length is not fixed, 
the body has a sharp nose. In the case of a body of given length and 
base diameter, the additional effect of curvature of the flow. over the 
.surface is investigated to determine its influence on the shape for min- 
imum drag. Th@^ .effect is. ‘ to increase the bluntness of the shape in the 
region of the hose and the curvature in the region downstream of the nose. 

Several bodies of revolution of fineness ratios 3 and 5, including 
the calculated shapes of minimum drag for given length and base diameter 
and for given base diameter and surface area', were tested in the Ames 
10— by 1^— inch supersonic wind tunnel at Mach numbers from 2.73 to 5-00. 

A comparison of theoretical and experimental foredrag coefficients indi- 
cates that the calculated minimum— drag bodies are reasonable approxima- 
tions to the correct shapes. It is verified, for example, that the body 
for a given length and base diameter has as much as, 18 percent less fore- 
drag than a cone of the same fineness ratio. 


INTRODUCTION 


i 


The shapes of nonlifting bodies of revolution having minimum pres- 
sure drag at supersonic speeds have been the subject of numerous theo- 
retical investigations. Karman (reference 1) determined the shape Qf 
such a body (neglecting base drag) with given length and base diameter. 
Somewhat later Haack (reference 2), Ferrari (reference 3 ) > Lighthill 
(reference 4), and Sears (reference 5) calculated body shapes having 
minimum pressure drag for various other given conditions using methods 
similar to those first employed by Karman. In all these investigations 
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the assumption of small perturbation, potential flow was made. It is 
to be expected, therefore, that the shapes obtained by these investi- 
gators are representative of minimum-drag body shapes of practical 
fineness ratios only at low supersonic Mach numbers . 

Perhaps the first calculation of the shape of a body having mini- 
mum drag was made by Newton (reference 6) using a method analogous to 
the present day calculus of variations. Newton was concerned with 
determining the body of given length and base diameter having minimum 
resistance when moving at sufficiently high speeds to insure that the 
inertia forces are large compared to the elastic forces in the immersing 
fluid. Thus, as shown by Sanger (reference 7 ) and Epstein (reference 8 ), 
the law of resistance adopted by Nevton approximates that (neglecting 
viscous forces) for hypersonic air flows. According to this law, the 
local resisting pressure is proportional to the square of the free— 
stream velocity component normal to the body surface. Legendre (see, 
e.g. , reference 9 ) further investigated Newton's problem and concluded 
that if no restrictions were imposed on the variation of slope along 
the surface, a body having a meridian curve composed of jagged lines 
(sharp edges forward) could be constructed which, with this law of 
resistance, would have less drag than Newton's body. It may easily be 
deduced, however, that Newton's law of resistance would not be satisfied 
on the surface of Legendre's body since gas would be trapped in a number 
of regions along the jagged contour. It may be shown in fact that when 
this law of resistance is satisfied at the surface — in which case the 
surface angles must lie between 0 and it/2 radians — then Newton's body 
may be considered the minimum pressure drag body for the given conditions. 

Aside from the investigations of Newton and Legendre, it appears 
that little, if any, work has been done which is applicable to the 
determination of body shapes which are especially suited, from the drag 
standpoint, to flight at high supersonic ‘ airspeeds. It has therefore 
been undertaken in the present report, using Newton's law of resistance 
and the calculus of variations, to determine body shapes having minimum 
pressure drag (neglecting base drag) at high supersonic speeds for 
various combinations of the conditions of given length, base diameter, 
surface area, and volume. In the case of a body of given length and 
base diameter, the effect of curvature of the flow over the surface is 
investigated to determine its influence on the shape for minimum drag. 

Several bodies of revolution, including two of the bodies deter- 
mined from this analysis, were tested at Mach numbers from 2.73 to 5.00 
in the Ames 10— by l4— inch supersonic wind tunnel. Foredrag data at 
zero lift obtained from these tests are compared with the analytic pre- 
dictions to assess the accuracy of the theoretical considerations. 
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local cross— sectional area of "body 
local speed of sound 
drag coefficient 
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pressure foredrag 
maximum "body diameter 
integrand function 
drag parameter! 
hypersonic similarity parameter f ) 

body length. 

Mach number 


(—) 

\2«q n / 
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distance measured normal to surface of body 

exponent in equation defining shapes of experimental test bodies 


pressure coefficient 

static pressure 
dynamic pressure 




o 


radius of curvature of streamline in plane containing axis of 
symmetry (i.e., meridian plane) of body 

body surface area 

resultant velocity 

body volume 


x,y coordinates of pointeoh meridian curve of body (origin of coordinate 
system coincides with nose of body, and x axis coincides with 
axis of symmetry) 


b 
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7 ratio of specific heat at constant pressure to specific heat at 
constant volume 

5 angle (in meridian plane) between free-stream direction and 
tangent to body s.urface 

X Lagrange multiplier 

p density 

Subscripts 

0 free-stream conditions 

1 values at nose point of meridian curve 

2 values at base point pf minimizing curve 

+ right-hand limiting value of quantity at corner on minimizing curve 
— left-hand limiting value of quantity at corner on minimizing curve 
B values along meridian curve 


THEORY 


The investigation undertaken here is concerned with the shapes of 
nonlifting bodies of revolution having minimum pressure foredrag at 
high supersonic airspeeds. Difficulties inherent in the calculation of 
these shapes make it desirable to simplify the drag equation insofar as 
is practicable, consistent with retaining the salient features of the 
dependence of drag on body shape and free-stream conditions. Likewise, 
in view of the several conditions to be treated (viz., given length, 
base diameter, volume, and surface area), it is convenient to set up a 
procedure of analysis to fit the general problem at hand. These fundar- 
mental considerations will be discussed prior to the determination of 
specific minimum-drag shapes^ 
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Fundamental .Considerations 


Simplified drag theory .— As pointed out in the introduction, 
Newton's law of resistance applies approximately to bodies traveling 
at high supersonic airspeeds. This observation has basis in the fact 
that at such speeds the inertia forces predominate over the elastic 
forces in the disturbed air. Thus, oblique shock flows approach the 
corpuscular- type flows treated by Newton as the Mach number of the free 
stream becomes large compared to 1. If it is further assumed that 
7 of the disturbed fluid approaches l, 1 the shock-wave angle approaches 
the flow-deflection angle (see sketch). In this case the pressure 


" 2 ^ 


coefficient at a point just downstream of the wave is given by the 
simple expression (reference 8) 

P = 2 sin 2 5 (1) 

This equation is recognized, of course, as being (aside from the con- 
stant multiplier) a mathematical statement of Newton's law of resist- 
ance for corpuscular or impact— type flow. 

When the curvature of the body, and hence of the disturbed flow, 
is small in the stream direction, equation (l) should also predict the 
pressure coefficients at the surface of a body since, in this case, the 
centrifugal forces in the thin layer of -air (sometimes referred to as 
the hypersonic boundary layer) between the shock and the surface should 
not appreciably alter the impact pressures. When the curvature of the 
body is large in the stream direction, centrifugal forces in the fluid 
between the shock and the surface may alter the pressures at the surface 

^Chis assumption was perhaps first suggested by Epstein (reference 8). 

It has the advantage of simplifying the impact force analysis without 
significantly influencing its accuracy. In other words, these forces 
in flows about axially symmetric bodies are relatively insensitive to 
changes in 7 from 1.4 to 1. Interestingly enough, present i nd ie ar- 
tions are that disturbed air flows in flight at extremely high Mach 
numbers will be characterized by values of 7 below 1.4 (see ref- 
erences 10 and 11). 
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from those just downstream of the shock. Busemann (reference 12) 
investigated this problem and found that the pressure coefficient at 
a point on the surface of a body' curved in the stream direction is 
given by the relation ^ 

P = 2 ^sin 2 5 + sin 6 ~ J' cos 5 dA ^ ( 2 ) 

' o 

in the limit as M— 5* »and 7 — > 1. 

In order to assess the accuracy with which the preceding equations 
may be expected to provide the pressure distributions, and thus pres- 
sure drags, on bodies operating at high supersonic airspeeds, the pre- 
dictions of these equations are compared in figure 1 with those of the 
method of characteristics (obtained from reference 13 ) for an ogive 
operating at a value of the hypersonic similarity parameter K (ratio 
of free-stream Mach number .to slenderness ratio) equal to 2, corre- 
sponding to a free-stream Mach number of 6 . It is evident that the 
theory of Busemann (equation (2)) yields far too low pressures down— , 
stream of the nose, while the simple impact theory (equation (l)) is in 
reasonably good over— all agreement with the method of characteristics . 2 
The relatively poor predictions of the Busemann theory are associated 
with the fact that it strongly overestimates centrifugal— force effects 
at free-stream Mach numbers which are large compared to 1, but for which 
7 of the air flow downstream of the bow shock is closer to 1.4 than 1 
(i.e., at flow conditions of principal interest in this paper). This 
matter will be discussed in greater detail later in the paper. Agree- 
ment comparable to that just discussed is obtained with the other 
results presented in reference 13 for K = 2. For lower values of K 
the agreement of the impact theory with the method of characteristics is 
somewhat poorer, as would be expected; however, it does not become 
unacceptably poor except for values of K below 1 (e.g. , the pressure' 
coefficients differ by from 0 to 35 percent for a K of l/2). It is 
therefore concluded that for values of K greater than 1 , equation (l) 
may be used with acceptable accuracy for the purposes of this paper to 
predict the pressure distributions and thus pressure drags on bodies . 3 

2 

The characteristics solutions of reference 13 were carried out for 
7 = 1.4. Free flight at the larger values of K considered in this 
report would produce values of 7 downstream of the bow shock slightly 
less than 1.4; however, the small decreases in 7 would not signifi— 

3 cantly alter the results presented. 

The approximate theory of reference 14 might also be employed in an 
analysis of this type — this theory is, in fact, more accurate than 
those discussed, so long as the flow is everywhere supersonic. How- 
ever, the theories considered here have the advantage, as will be 
apparent later, of predicting approximately correct values of surface 
pressures for arbitrarily large surface slopes. 
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For this reason, and "because of its simplicity, it is employed through- 
out the subsequent analysis. 

If the manner in which the pressure coefficient varies over the 
surface is known, it is a simple matter, of course, to evaluate the 
pressure drag of a body. Neglecting the base-drag contribution, we 
have then 

CL q ltd 2 p l 

D = -2 — = 2* q Q J Pyy’dx (3) 


where y' denotes the derivative dy/dx, This equation may be 
expressed in a form more convenient for use here 


I 3) = 


D 


2*q r 


Pyy'dx 


(4) 


If P in this expression is replaced by its value given in equation (l) 


(noting that sin 2 6 = 


_ y 


l2 


l+y' 


= 


r ), there is then obtained the relation 


l 3 

2yy’ 


dx 


l+y* 


(5) 


It remains now to consider the procedure for employing this expression 
in combination with the methods of the variational^ calculus in order to 
determine the desired minimum-drag body shapes. 


Procedure for Calculating Minimum-Drag Bodies 


The calculation of minimum-drag body shapes of interest here is 
equivalent to determining the form of the function y = y(x) which 
minimizes the integral defined in equation (5) for the various given 
conditions. In considering the procedure for carrying out this calcu- 
lation, however, it is convenient, for reasons that will be apparent 
later, to Write equation (5) in a form which effectively yields the 
total drag as the sum of the drag on any finite region of infinite slope • 
at the nose plus the drag on the surface downstream of the nose. Thus 

we have ' „ 

x 2 


I D = 7i 2 + 



SZZlldx 

l+y* 2 


( 6 ) 
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where th6 variable limit x 2 is introduced to permit variations in 
body length. The conditions of given volume or given surface area are 
fixed "by the auxiliary requirements that, respectively, 

y 2 dx = const. (7) 


x 2 



or (neglecting base area) 




dx = const. 


( 8 ) 


When the length and hajse diameter are given, the problem is simply 
to minimize the function Ip given by equation (6). However, according 
to the isoperimetric rule of the calculus of variations (see, e.g., 
reference 15), the problem of minimizing the function subject to 

the auxiliary condition given by equation (7) or (8), is equivalent to 
minimizing the new function J D > where 

j D - h + x i (? > 

or 

J D = I D + ^ ( 10 ) 

depending on whether the volume or surface area is given. The parameter 
X is a constant, sometimes called the Lagrange multiplier. 

With the aid of equations (6) through (10), the integrand functions 
to be minimized can be immediately written. These 'functions are as 
follows : 

case a, given length and base diameter 

f = 2 yy ' 3 
l+y* 

case b, given volume and length or base diameter 

f = gXXl! + ^ 
l+y' 2 


( 11 ) 


.(12) 
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case c, given surface area and length or base diameter , 

f = g yy .— + Xyjl + y ,a ( 13 ) 

l+y' 

Now any function y = y(x) which minimizes equation ( 6 ), ( 9 )> or 
(10)' must, irrespective of the given conditions, satisfy the Euler 
equation (for zero first variation of Ip or Jp with small changes in 
the function y(x)) 

■'sV- f y-o < l4 > 

where f v ' and f v denote the partial derivatives and respec— 

j J 67 dy 

tively. Since the integrand functions given above are free of the 
independent' variable, the first integral of the Euler equation for 
these functions follows immediately, namely, 

y* fyt — f = const. ( 15 ) 

Substituting, successively, equations (ll), (12), and (13) into this 
equation there are then obtained the expressions 


^yy ,a 

(l+y ,2 )‘ 


= const. 


and 


%y' 

(l+y' 2 ) 


— Xy 2 = const. 


y 3 ; _x 

(l+y * 2 ) 2 */l+y' 2 - 


= const. 


(16) 

(IT) 


( 18 ) 


for cases a, b, and c, respectively. Solutions to these differential 
equations satisfying the terminal conditions on the bodies are minimiz- 
ing curves for the given conditions. 

When the end points of a minimizing curve are not fixed, other 
terminal conditions must be imposed on the function y = y(x). Thus, 
to determine the ordinate at the nose, it is required that (see ref- 
erence l 6 ) 


C 
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L 7» dy 72 ) 


for cases a and 1) , while 


= 2y 


y=yx 


(i+y’ 2 ) 


0 


( 19 ) 


y=yi 


L f y - ^ (y2 * x f J 


= yi 


y=yi 


2 iL±2 Ill + -£L-(2 + \) 


L (l+y* 2 ) Vl+y'‘ 


= 0 


-r=yi 


( 20 ) 


for case c. Similarly, when the length is not given it is necessary 
that 


y’fy. -f 


= 0 


( 21 ) 


X=X 2 

and vhen the hase diameter is not given it is required that 


f t 

7 


x=x 2 


= 0 


( 22 ) 


In addition to the above described conditions, two checks must be 
made to determine completely the shape of a minimizing curve. The first 
of these checks entails ascertaining whether there are any corners 
(between the end points) on the curve. This is accomplished by deter- 
mining whether the function y = y(x) can satisfy the requirement that 
(see reference 15 ) 

f = f 

y*+ y* 

at a point of discontinuity in y*. If this equation is not satisfied, 
no corners exist. The second check requires that the Legendre condition 
(for a positive second variation), 

Vy - 0 (2k) 

be satisfied everywhere on the curve. With the aid of these checks, the 
minimizing curves for various combinations of the conditions of given 
length, base diameter, volume, and surface area can be uniquely defined. 
The calculation of these curves for several such combinations is now 
undertaken. 


( 23 ) 
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Calculation of Minimum-Drag Bodies 


Given length, and base diameter .— Equations (l6) and (19) give the 
first integral to Euler’s equation and the terminal condition at the 
nose, respectively, for these given conditions. It is evident upon 
examining these equations that the minimizing curve cannot, in general, 
pass through hoth the points (o,o) and (x 2 ,y 2 ), "but must, in fact, have 
its forward termination point at (o,yi.) with y^ = 1. With this infor- 
mation, the minimizing curve 
namely, 

y 


x 


can he represented in parametric form, 

.2 


.7i (i + y* ) 

” 5 

4 y» 


3 + ^_l + zny. 

4 


yi t 

4 ^4y’ 4 y’ 2 


(25) 


It is easily shown with the solution to the Euler equation and equa- 
tion (23) that there are no corners on the minimizing curve; 4 thus the 
variation of y with x is readily determined with the relations of 
equation (25) for a given 2 and d (corresponding to a given x 2 and y 2 ) 
of a body. These relations for a body of given fineness ratio can be 
shown to be equivalent to those originally developed by Newton (see ref- 
erence 6) . 


Given length and volume .— For these given conditions, the terminal 
conditions (equations (19) and (22)) require the slopes at the nose and 
at the base to be, respectively, yi’ = 1 and y 2 ’ =0. The first inte- 
gral to the Euler expression (equation (17)) then leads to the following 
parametric representation of the minimizing curve: 


3 

2y» 

X(l+y’ 2 ) 2 



7i - kyi 2 

x 





( 26 ) 


Similarly, it can be shown that there are no comers between (o,yi) 
and (x 2 ,y 2 ) on any of the minimizing curves to be treated here. 
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From the relations of equation (26) it is clear, again, that the mini- 
mizing curve cannot pass through (o,o), the condition yi* = 1 deter- 
mining a value yi >0. These relations, together with the volume 
condition (equation (7)) and the given length condition, serve to 
determine yi and A and thus, of course, the shape of the entire "body. 
As the length approaches 0, A becomes infinitely negative; while, as 
the length becomes infinitely large, A approaches 0. (in the latter 
case the body shape approaches the minimum-drag shape for the given 
length and diameter condition, Z/d -$>».) Intermediate negative values 
of A correspond to intermediate values of length for. a given volume. 


Given length and surface area .— In this case a first integral to 
the Euler equation is given by equation (l8), and the parametric repre- 
sentation of the minimizing curve may be written immediately in the 
form 


const, (l+y* 2 ) 

7 3 2 3 / 2 

4y* — A(l+y' ) 


> 



(27) 


Upon examination of this equation and equations (20) and (22), it 
becomes apparent that, again, the minimizing curve cannot go through 
the point (o,o). The latter equations determine uniquely, however, the 
values of yx’ (yi’< 1), and y 2 * (o<y2*<yi*)> in terms of the param- 
eter A. Similarly, the length and surface— area condition in combina- 
tion with the above equations determines the value of A. Thus it is 
easily shown that the practical range of A is from -2 to 0 (corre- 
sponding to body lengths of from zero to infinity for a given surface 
area — in the latter case the Newton body is again obtained),' 


Given base diameter, and volume .— With these given conditions, the 
first integral to the Euler relation is given by equation (17) , while 
the terminal conditions at the fore-and-aft ends of the body are fixed 
by equations (19) and (21), respectively. It is evident that the mini- 
mizing curve must, in general, pass through the origin in order to 
satisfy all these equations in addition to the Legendre condition 
(equation (2k)). The shape of the minimizing curve may thus be defined 
parametrically as follows: 
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X (l+y ,2 )‘ 


2 y* + 3y* ' 
^ (l+y‘ 2 ) 2 


where yj.' = 0. Combining this expression with equation (7), there is 
•then obtained for the volume of the body 


V = (y 2 ' 4 + 6y 2 ' 2 + ^5) 

120 y 2 ' 


The range of X for which these results are applicable is from 

zero to 3 «/3/4y 2 , corresponding to a volume range from infinity to 

rt yp 3 -s/3/5. For a given y 2 and a given V>rf y 2 3 */375 (corresponding 

to ^> 1 / 3 / 2 ), equation (29) has two solutions in y 2 '. One solution 

yields values of y 2 ' greater than >J 3 , which result violates the 
Legendre condition (see equation (24)), while the other yields permis- 
sible values less than »/3. When y 2 and y 2 ' are known, X may then be 
determined from the first relation of equation (28), namely, 

X. =- 72 ' 3 5 (30) 

72 (i+yp* 2 ) 

The determination of y and x follows directly, of course, from 
equation (28). 

Given base diameter and surface area .— In this case equations (18), 
(20), and (21) determine the shape of the minimizing curve as being 
simply a straight line 

/ (xA , 2/3 , , 

y = x / : £75- (31) 

J 1 - (xAf /3 


where the parameter X is given by the equation 


5 The solution given here is not applicable to bodies of extremely small 

fineness ratios (viz., as can be easily deduced from 

equation (28). 
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X = 4 (rty 2 7s) 3 (32) 

Thus, the mininnmt-drag body for given base diameter and surface area is 
a cone. 


Comparison of Minimum-Drag Body Shapes 


The previous calculation of minimum-drag bodies reveals, two 
general characteristics of their shapes; namely, when the length is 
given (fixed) the bodies assume blunt noses, whereas, when the length 
is not given (i.e,, is free), the bodies assume sharp noses. The 
former characteristic may be traced to the fact that with the length 
restricted, the net drag is reduced by accepting higher pressures on a 
relatively small area of large slope near the nose, thus achieving 
lower pressures on a relatively large area of small slope near the base. 
On the other hand, when the length is not restricted it is evident that 
a sharp rather than a blunt nose will obtain for minimum drag, since 
the drag of any blunt-nosed body can be reduced by simply relaxing the 
requirement on length, thereby allowing the body to be made sharp nosed 
and generally more slender. 

In order to permit a quantitative comparison of the shapes of the 
calculated minimum^drag bodies, typical meridian curves for these bodies ' 
are shown in figure 2. For simplicity the bodies are compared on the 
basis of the same fineness ratio — ordinates have been plotted to an 
expanded scale to better indicate the relative shapes. The maximum 
bluntness is evidently obtained when the drag is minimized for a given 
length and surface area, while the maximum sharpness (a cusp nose) is 
obtained when the base diameter and volume are given. It is apparent 
from figure 2 that the flat-nosed portions of the meridian curves for 
the given length bodies are in all cases very small. For example, 
yi equals 0.0050y e for the body of given length and volume. On the basis 
of several calculations it is indicated, as might be expected, that the 
degree of bluntness will increase with decreasing fineness ratio. 

It is also of interest to compare minimum-drag body shapes deter- 
mined with the aid of the linear theory (see, e.g., reference 2) with 
those found using the impact theory, that is, bodies especially suited 
for flight at low and high supersonic speeds, respectively. Such a 
comparison is shown in figure 3 for the case of given length and base 
diameter. It is seen that qualitatively the shapes are similar although 
the minimum-drag body for low supersonic speeds is generally the fatter 
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of the two. 6 Comparisons of the results of this paper with those of 
reference 2 for other given conditions also indicate qualititive agree- 
ment as to general "body shapes despite the marked difference in the 
laws governing the surface pressures. 

All the preceding analysis has "been predicated on the assumption 
that the flow of air at high supersonic speeds may, insofar as pressure 
forces are concerned, be approximated by a Newtonian- type flow. It 
remains now to test the accuracy of this assumption and other aspects 
of the analysis by experiment. 


EXPERIMENT 


It has been undertaken to obtain a partial check on the findings 
of the preceding theoretical analysis by determining experimentally the 
foredrags on a family of bodies, including calculated minimumr-drag 
shapes, at Mach numbers from 2.7 to 5.0. The analysis may be expected 
to apply at least approximately in this range since the corresponding 
valuesof the hypersonic similarity parameter K were, for the most 
part, greater than 1. A brief description of these tests is presented. 


Test Apparatus 


Wind tunnel .— The tests were conducted in the Ames 10— by 1 4— inch 
supersonic wind tunnel which is a closed test section, nonreturn tunnel 
of the continuous— flow type. A schematic view of the wind tunnel and 
auxiliary equipment is shown in figure 4. The 10— inch dimension cor- 
responds to the width of the tunnel which is constant from Just upstream 
of the nozzle to downstream of the diffuser. The 1 4— inch dimension is 
the nominal value of the depth of the test section. The top and bottom 
walls of the wind tunnel are rigid (as are the sidewalls) and form the 
nozzle and the converging-diverging diffuser. These walls are trans- 
lated and/or rotated relative to each other during operation of the 
tunnel in order to vary the test— section Mach number. At the same time 
the diffuser shape is adjusted to decrease the intensity of the terminal 
shock system, thus decreasing the power required to operate the t unn el. 

^Part of this .difference in shapes may stem from the fact that the body 
derived using linear theory was required to have zero slope at the 
base. Also, as will be shown later, the true minimum-drag shape at 
high supersonic airspeeds may be somewhat fatter than that obtained 
using inpact theory, due to the fact that centrifugal forces are neg- 
lected in this theory. 
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At present the tunnel' is operated at test Mach numbers from 2.7 to 5.0; 
operation at higher Mach numbers is prohibited by excessive, condensation 
of the air in the test section. 

Flow in the diverging portion of the diffuser is stabilized by 
removing part of the low-energy boundary— layer air through scoops 
located at the minimum section of the diffuser. This air is diverted 
to the vacuum pumps shown beneath the wind tunnel in figure 4. The 
remainder of the air is exhausted either directly to the atmosphere, or, 
when higher pressure ratios are required, through the air-driven, two- 
stage, centrifugal-evacuator system shown at the lower right of figure 4. 

Supply air for the wind tunnel and. drive air for the centrifugal 
evacuators are provided by two stages of centrifugal compressors (not 
shown in fig. 4) at absolute pressures up to a maximum of 6 atmospheres. 
The supply air is passed through coolers and silica— gel dryers to reduce 
the absolute humidity to less than 0.0001 pound of water per pound of 
dry air. The drive air is piped to the evacuators from immediately 
downstream of the second stage compressors. 

The schlieren apparatus used for visual investigations of the flow 
is shown schematically in figure 4. The system consists basically of a 
light source, a spherically ground mirror with a radius of curvature of 
50 feet, an adjustable knife edge, and a ground glass viewing screen 
which may be replaced by photographic film packs. Two plane mirrors 
are used to reduce the length of the system while retaining the 50— foot 
light path. The entire apparatus is installed in a light-tight enclo- 
sure straddling the tunnel test section. The test section is traversed 
twice by the light beam, thus increasing the sensitivity over that of 
the more conventional types of schlieren systems. 

The schlieren apparatus discussed previously was converted to a 
shadowgraph system for the tests reported here in order to obtain 
clearer outlines of the bow shock waves which were of primary interest. 
This conversion was accomplished by simply replacing the spherical, 
mirror in the schlieren system with photographic plates and operating 
the light source as a spark — in this case the test section is trav- 
ersed only once, of course by the light beam. 

In figure 5 are shown the variations in Mach number along the test- 
section center line as determined from pitot and static pressure meas- 
urements made at 1.5— inch intervals in the axial direction. These var— - 
iations do not exceed +0.02 except at Mach numbers of 4 and 5 where they 
amount to ±0.03 and ±0.04, respectively. The corresponding maximum 
variations in- pressure coefficient are ±0.005, 0.003, and 0.002, res- 
pectively. 
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The variation of Reynolds number per foot with Mach number for the 
present tests is shown in figure 6. 

Ins t rument at i on . — Aerodynamic drag was measured by a single- 
component strain— gage balance. A .cutaway drawing of the balance is 
shown in figure 7. Tare forces on the sting supporting the model were 
eliminated by enclosing the sting within a shroud that extended to 
within 0.o4o inch of the model base. 

Pressures at the model base were determined with the aid of a low- 
pressure McLeod gage. Stagnation pressures in the reservoir were 
measured on a Bourdon— tube gage. 

Models .— Five models of fineness ratio 3> and three models of 
fineness ratio 5 were tested. With the exception of an Z/ d = 3 tan- 
gent ogive (this shape was included as being typical of those in common 
usage), all models had meridian section shapes given by the equation 



where n was given values of 1, 3/4, l/2, and l/4. When n = 3/4, 
the body shapes defined by the above expression closely approximate the 
minimum-drag shapes for given length and. base diameter (equation ( 25 )) 7 
for z/d = 3 and 5 (see fig. 8). 

When n = 1, the cone is, of course, obtained which is the minimum- 
drag body for a given base diameter and surface area. Minimum-drag 
shapes for two different given conditions are thus included among the 
bodies tested. 

Photographs of the eight models tested are shown in figure 9. The 
Z/d = 3 bodies (fig. 9(a)) are, from left to right in the photograph, 
the cone, 3/4-power body, l/2-power (parabolic) body, l/4— power body, 
and the tangent ogive which has a profile section radius of curvature 
of 9.25 body diameters. From left to right in figure 9(b) are the 
Z/d = 5 cone, 3/4— power body, and l/2-power body. The base diameter 
of all models was 1 inch. 


The accuracy of this approximation increases with increasing values of 
Z/d as can easily be seen upon examination of equation ( 25 ). 
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Accuracy of Test Kesults 


The accuracy of the foredrag coefficients is effected by uncer- 
tainties in the measurements of the following quantities: stagnation 

pressures, free-stream pitot pressures, free-stream static pressures, 
base pressures, and the forces on the models as measured by the strain- 
gage balance. 

Both static and free-stream dynamic pressures were determined from 
wind-tunnel calibration data and stagnation— pressure readings. The 
latter measurements were accurate to within ±l/2 percent, thus static 
and dynamic pressures are uncertain by this amount plus possible cali- 
bration errors of ±1 percent over the Mach number range of -the tests. 

The uncertainty in foredrags due to inaccuracies in the determination 
of base pressures does not exceed ±1 percent. 

Because of the small drag forces measured, the source of greatest 
error was the strain-gage balance system. The uncertainty in drag due 
to zero shifts, thermal effects, and friction varied from approximately 
±2 percent at the lower Mach numbers to ±6 percent at the highest Mach 
numbers . 

The models were tested on the tunnel center line between stations 
46 and 54. In this region the variations in stream static pressure 
coefficient discussed previously were too small to warrant the applica- 
tion of a bouyancy correction to the measured drags. 

The combined effects of all the sources of error result in probable 
uncertainties in measured foredrag coefficients of from ±0.001 at tbs 
low Mach numbers to ±0.005 at a Mach number of 5.O. 8 

In order to reduce this error in the data presented here, partic- 
ularly at the higher Mach numbers, several measurements were made at 
each Mach number and the average values of foredrag coefficients were 
employed. 


RESULTS AND DISCUSSION 


The - variations with Mach number of the pressure drag and skin- 
friction drag for cones with laminar boundary layers can be accurately 

Q F — _ --- _ - 1 --- - 

At Mach numbers near 5 there is some additional uncertainty due to the 
presence of a small amount of condensed vapor (including O 2 and N 2 ) in 
the test— section air stream; however, as will be evident later, this 
uncertainty is inconsequential to the present tests. 
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predicted with available theories. The Eeynolds numbers of the present 
tests were sufficiently low to insure laminar boundary layers in all 
cases (see fig. 10 showing sharply defined boundary layers in -shadow— 
graphs of the flow about the l/d = 3 cone and l/2— power body). .Since 
such a calculation provides a reliable check on the accuracy of the 
experimental results of the present tests, it has been carried out for 
the cone of l/d = 3. The pressure drag was determined from ref- 
erence 17. The skin— friction drag for a flat plate was determined from 
the results of reference 18 and modified by the transformation of 
Hantzche and Wendt (reference 19) for the case of the cone. A compari- 
son of the foredrag coefficients for this cone as determined experi- 
mentally and by the calculations described above is shown in figure 11. 
Good agreement between experiment and this relatively exact theory is 
observed over the Mach number range of the tests, thus substantiating 
the accuracy of the experiments. The increase in foredrag coeffi- 
cients at Mach numbers above 4.5 is due to the increase in friction 
drag which in turn results from the rapid decrease in Eeynolds number 
with increasing Mach number in this range (see fig. 6). 

The variations with Mach number of the foredrag coefficients of 
the eight test bodies are shown in figure 12. With the exception of 
the l/4— power body, these variations are similar to those for the fine- 
ness ratio 3 cone previously discussed. In the case of the l/4— power 
body Cp increases with increasing Mach number over the entire test 
range. This behavior may result from the fact that, while the skin- 
friction drag coefficient is increasing with Mach number, the pressure 
drag coefficient for this very blunt body, in contrast to the more 
slender bodies investigated, probably decreases little if at all. 9 

The results presented in figure 12 show that, as predicted in the 
analysis, the 3/4— power bodies have the minimum foredrags of all bodies 
tested having the same fineness ratios. Percentagewise the improve- 
ment is larger at the higher Mach numbers. This result might be antici- 
pated from the simplified drag theory used in the analysis (this theory 
is strictly applicable only at high Mach numbers for which the value of 
the hypersonic similarity parameter K for the flow is greater than 1). 
Except in the case of the bluntest body, however, the simplified theory 
yields foredrags that are appreciably lower than those obtained experi- 
mentally at all Mach numbers. Nor would the addition of the correct 
skin— friction drag to the theoretical pressure drag bring theory and 
experiment into agreement since the theoretical drag is generally too 
low. The differences between the experimental and theoretical results 
are illustrated in the following table where the drag coefficients are 
compared for all the bodies at a free-stream Mach number of 4.0. 

®In the limit as the body becomes completely blunt, that is, a cylinder, 
the pressure drag coefficient increases slightly with increasing Mach 
number in this range. 
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Model 

Eoredrag coefficients 

Experiment 

Equation 1 

l/d = 3 first power (cone) 

0.079 

0.0541 

Z/d = 3 3 A power 

.067 

.0455 

z/d = 3 l /2 power 

.085 

.0692 

Z/d = 3 lA power 

.249 

.2878 

Z/d = 3 ogive 

.096 

.0705 

z/d = 5 first power, (cone) 

.043 

.0198 

Z/d = 5 3A power. 

.037 

.OI 67 

Z/d = 5 l/2 power 

.046 

.0300 


At higher Reynolds numbers and higher Mach numbers corresponding to 
larger values of K (the highest value of K in the present tests was 
I. 67 ) it would he expected that theoretical and experimental foredrag 
coefficients would he in better agreement. 

A check on the overfall accuracy with which the optimum shapes are 
predicted by the analysis is obtained by comparing theoretical and 
experimental values of the relative foredrag coefficients of the test 
bodies. Such a comparison is given in figure 13 where the ratios of 
the foredrag coefficients of a test body to the corresponding coeffi- 
cients of the cone of the same fineness ratio are shown as a function 
of the exponent n in equation ( 33 ) defining the shapes of the test 
bodies. The theoretical predictions of the impact theory appear to be 
in good agreement with the experimental results at the higher values of 
n (approximately n>0.6). Thus it is suggested that the 3 /4-power body 
is a reasonable approximation to the correct minimum foredrag shape of 
given fineness ratio. At the lower values of n, however, it is indi- 
cated that in contrast to absolute drag, previously discussed, the 
relative drag is significantly overestimated by this theory. This 
result is not entirely surprising since the theory neglects centrifugal- 
force effects in the disturbed flow, and these effects must appreciably 
alter the pressures over the highly curved noses of the blunter bodies 
(see fig. 9 ). 

As discussed earlier, the Busemann theory for infinitely high Mach 
numbers overestimates these effects at the Mach numbers of interest here. 
It has therefore been undertaken in appendix A of this paper to obtain 
a better estimate of centrifugal forces by accounting approximately for 
the decrease in these forces (at finite but high Mach numbers) associ- 
ated with the increase in the lateral extent of the disturbed flow field 
with increasing distance downstream from the nose of the body (see 
fig. 10). The predictions of the modified impact theory shown in 
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figure 13 were obtained with the aid of this estimated centrifugal-force 
effect (see equation (A9)) in combination with equations (l) and (3). 

It is indicated that this theory is' markedly superior to the impact 
theory at the lower values of n, corresponding to the blunter bodies, 
over the test Mach number range. The estimate of the centrifugal 
forces would thus appear. to be in fair agreement with the actual magni- 
tude of these forces. 

It is also indicated in figure 13 that at the highest test Mach 
numbers (hence highest values of K) the modified theory is generally 
somewhat superior to the impact theory. This result suggests that 
improved approximations to the correct minimum-foredrag shapes for 
values of K appreciably greater than 1 may be obtained by using this 
theory rather than the simple impact theory. As a first step in the 
direction towards obtaining these shapes, the minimum-drag body of 
given length and base diameter has been calculated with the aid of the 
modified impact theory. This calculation is carried out in appendix B, 
and the resulting shape for a given fineness ratio is shown in figure l4. 
Newton’s body of the same fineness ratio is also shown for comparative 
purposes and, as would be expected, this body is somewhat less blunt 
than the former shape in the region of the nose and has less curvature 
in the region downstream of the nose since centrifugal relieving effects 
on the surface pressures are not considered by the impact theory. 


CONCLUDING REMARKS 


It has been undertaken in this report to determine approximately 
the shapes of several bodies having minimum pressure foredrag at high 
supersonic airspeeds. With the aid of Newton’s law of resistance and 
the calculus of variations, an investigation was carried out for various 
combinations of the conditions of. given body length; base diameter, sur- 
face area, and volume. In general, it was found that when the length is 
fixed, the body has a blunt nose (i.e., a finite area of infinite slope 
at the nose) as in the classical problem considered by Newton; whereas 
when the length is not fixed the body has a sharp nose. 

Several bodies of revolution of fineness ratios 3 and 5> including 
the calculated minimum-drag bodies for given length and base diameter 
and for given base diameter and surface area, were tested at Mach 
numbers from 2.73 to 5.00 in the Ames 10 — by 1 4— inch supersonic wind 
tunnel. A comparison of the relative theoretical and experimental fore- 
drag coefficients indicated that the calculated minimum-drag bodies were 
reasonable approximations to the correct shapes. It was verified, for 
example, that the body for a given length and base diameter has as much 
as 18 percent less foredrag than a cone of the same fineness ratio. 
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The cone is, however, the calculated minimum-drag body for a given base 
diameter and surface area. 

The comparison between theory and experiment also indicated that 
the centrifugal forces in the flow about bodies curved in the . stream 
direction may significantly influence their drag. The relative extent 
of this influence was found to be predictable, particularly at the 
higher Mach numbers, with a simple modification to the impact theory of 
Newton. It was therefore suggested that improved approximations to 
minimum foredrag shapes at high supersonic airspeeds (for which the 
hypersonic similarity parameter has a value appreciably greater than 1) 
may be calculated with the aid of the modified impact theory. Such a 
calculation was carried out for the body of given length and base diam- 
eter. The resulting shape was found to be somewhat blunter in the 
region of the nose, and to have more curvature in the region downstream 
of the nose than the corresponding shape obtained using the simple 
impact theory. 


Ames Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Moffett Field, Calif. 
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APPENDIX A 

ESTIMATED EFFECT OF CENTRIFUGAL FORCES 
ON SURFACE PRESSURE COEFFICIENTS 


An estimate of the effect of centrifugal forces on the pressures 
at the surface of a hody operating at high hut finite Mach numbers may 
he obtained hy comparing the disturbance flow fields at these Mach 
numbers with that associated with infinitely large Mach number. 

At high Mach numbers the disturbed air flows in a relatively thin 
region (sometimes termed the hypersonic boundary layer) between the bow 
shock wave and the surface of the body (see sketch). 



The change in pressure from the surface to the shock due to centrifugal 
forces in the fluid is given by the equation 

Ap = / ^ dN = f dN 

' Jo- W Jo K 

assuming the directions of the normals to the streamlines between the 
surface and the shock do not differ appreciably from the direction of 
the normal to the surface. This expression is more conveniently written 
in the form 

- r*' 

Ap = 2 I p UdN (Al) 

R J o • 

where U and R are mean values of the velocity and radius, respec- 
tively, in the interval N. Now the mass m of air between the surface 
and the shock flowing (in unit time) by a point on the body is given by 
the relation 

r * 

m = 2 Tty / p UdN = * y 2 p Q U Q (A2) 


24 


NACA EM A51E27 


Combining equations (Al) and (A2) there is then obtained for the pres- 
sure change 

Ap = S | P 0 U 0 
R d 

or in coefficient form 

AP 


Now in the limit as the Mach number approaches infinity and y of the 
disturbed fluid approaches 1, the thickness of the layer becomes infini- 
tesimal and hence 

. ® = r B (A4) 

Similarly, it is easily shown (e.g. /with the compatibility equations 
applying along characteristic lines in axially symmetric supersonic 
flow) that - 

dU = 0 


Ap _ y U 
q ° R U 0 


(A3) 


along any streamline downstream of the bow shock, and thus that 

rJ 

ft 2U ° P 

u = — g- / y COS 5 dy 

y J o 


(A5) 


Hence, in this limiting case, equation (A3) takes on a form equivalent 
to that first deduced by Busemann (see second term on right of equa- 
tion (2)), and later derived in reference 20, namely, 


AP = 



y 

y cos 5 dy 


(A 6) 


where 


-L- = sin 6 ^ 

^ dy 


On the other hand when the Mach number is finite, but high, and -y 
of the disturbed fluid is closer to 1.4 than 1, the preceding evalua- 
tions of R and U are in considerable error since the hypersonic 
boundary layer, although thin, is no longer of infinitesimal thickness. 
This change in the boundary layer results from the fact that the bow 
shock is detached (except perhaps at the nose) from the surface of the 
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■body, the lateral distance from the surface to the shock increasing 
with increasing distance downstream from the nose (see sketch). Thus, 
for example, E would he expected to approach Eg only near the nose, 
while with increasing distance downstream of the nose it would he 
expected to become larger than Eg. From the pressure distributions 
presented in reference 13 it is indicated, _in fact, that for K>1 (the 
range of K’s of interest in this paper) E>>Bb near the maximum 
ordinate of the body. (This indication follows from the small values 
of the pressure coefficients near the maximum ordinate. ) ‘It is sug- 
gested, therefore, that at the high supersonic speeds under consider- 
ation, an approximation to E is given by the relation 



(AT) 


Similarly, in the case of TJ it no longer follows that the magni- 
tude of the velocity must be constant along streamlines downstream of 
the bow shock since pressure disturbances can now be_transmitted across 
streamlines. Thus a better first approximation to U than that given 
by equation (A5). may be obtained from the simple corpuscular or impact 
theory, 10 namely. 


U = U 0 cos 5 


(A8) 


Combining equations (A3), (A7), and (A8), the estimated change in pres- 
sure coefficient at the surface of a body due to centrifugal forces in 
high supersonic speed flow is obtained in the form 


or 


AP = -2- (1 _ JL) cos 5 

% 7a 

AP = £ (i _ L.) (sin 2 5) (A9) 

2 72 dy 


°With this theory, acceleration of the flow along a streamline is 
qualitatively accounted for. 
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APPENDIX B 

CALCULATION OF MINIMUM-DRAG BODY OF GIVEN LENGTH AND 
BASE DIAMETEB, WITH CONSIDERATION OF CENTRIFUGAL 
FORCES IN THE DISTURBED FLOW FIELD - 


For the purpose of this calculation, equations (l) and (A9) for 
the pressure coefficient are combined with equation (4) to yield the 
drag parameter in the form 


In = 7i +T (7i) + 


2 sin 2 6 + 2. (1_ X.) X s in 2 5 yy'dx (Bl) 


The term Jx represents the drag on any finite region of infinite 
slope at. the nose, while the function cp(yi), given by 11 


<p(yi.) - - (1 - ||) cos 2 Si 

2 72 


represents a "leading-edge thrust" due to the acceleration of the air 
flow about a comer (if it exists) at (o,yi). 

The expression (Bl) may be put in the form 

Id = yi 2+ / \'L 2 sin2 6 + 2 h sin2 5 j yyf “ dx 


whereupon (recalling sin 2 5 = — — ~) the integrand simplifies to a 

1+y* 

function f given by the relation 


1 This function may be obtained by evaluating 


yi+ e 


lim 
e-> 0 


yi- e 


T (1 - & ay (s “ 8)4y 


along the body- surface streamline about the corner at (o,yi)* 
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f = 77' 2 - 



With, the aid of this expression and equations (15) and (19) the para- 
metric representation of the minimizing curve can "be obtained in the 
form 



where 


and in general 


(B2) 


7i* = 1 


7i ± 0 


The minimizing curve given by these relations, similar to the curve 
obtained from the impact pressure treatment, does not have a comer 
between the points (o,yi) and (x 2 , Jz)* The minimum-drag shape defined 
by equation (B2) is compared in figure l4 with that determined earlier 
by considering impact pressures only. 
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Percent Body Length 



Figure tr Comparison of approximate and exact pressure distributions 
over a tangent ogive of fineness ratio 3 operating at a Mach 
number of 6 (K=2). 





Figure 2.— Minimum drag bodies for various given conditions. (t/d=5.0 for all bodies) 
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Figure 3.-Comparison of minimum drag bodies of given length and base diameter determined 

by linear theory and by impact theory. 
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Figure 4- Schematic view of the !0- by 14-inch Supersonic Wind Tunnel. 
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Figure 5.— Mach number distributions on the center line of the iO- by 14- inch Supersonic 

Wind Tunnel. 
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Figure 7. - Cutaway view of the drag balance used in the present tests. 
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Figure 8. - Comparison of profiles of minimum drag bodies of revolu- 
tion for given lengths and base diameters with the approximate 
profiles employed in the present tests. 




(a) Fineness ratio 3 bodies. 

Figure 9*- Photographs of the seven test bodies, the shapes of which are given by the equation 

y = &■ (x/l) n , and the l / d = 3 tangent ogive. 
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n= I (cone) 



(b) Fineness ratio 5 bodies. 


Figure Concluded. 
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M 0 = 4.01 


! 



(a) 2/d = 3 parabola M Q = 5-00 (b) 2 /d = 3 cone 


Figure 10.- Shadowgraphs of fineness ratio 3 conical and parabolic bodies at 
Mach numbers 2.73; 4.01, and 5*00 in the Ames 10- by l4-inch Supersonic 
Wind Tunnel. 
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Figure //. —The variation with Mach number of the foredrag coefficient at zero lift 

of the l/d s 3 cone. 
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(b)l/d= 5 bodies 

Figure i3r The ratios of fore drag coefficients of test bodies to 
foredrag coefficients of cones as functions of the exponent, n, 
in the equation defining body shapes. 
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Figure 14- The effect of centrifugal forces on the shape of the minimum drag body of given length 

and base diameter. (l/d°6.i8) 
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